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.^ ■ In this paper, a study of two-photon and two-gluon decays in the context of p-wave heavy 

^^ . quarkonia is presented. Within the covariant light-front framework, the annihilation rates of scalar 

and tensor quarkonium states are derived. In the absence of free parameters in this case, the 
Q^' results for the charmonium decay widths are consistent with the experimental data. However, in 
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comparison to other theoretical calculations, there are large discrepancies in our results regarding 
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I. INTRODUCTION 



Heavy quarkonium provides for a unique laboratory to study quantum chromo dynamics 
(QCD) regarding the bound states of a heavy quark- ant iquark system. Notably, the two- 
photon and two-gluon annihilation rates of p-wave heavy quarkonium are helpful for better 
understanding the details of quark-antiquark interaction and can function as stringent tests 
for a potential model. Regarding experimentation, the two-photon decay width of Xcj has 
been measured by many laboratories [l|, |2] and a new CLEO measurement was reported 



recently [4I . Regarding theory, relevant decay rates were first obtained through nonrelativis- 
tic approximation [31,14]; the relativistic corrections were included within the Bethe-Salpeter 



equation [5|-l8|, potential model |9|, relativistic quark model |l0l-ll2|. nonrelativistic QCD 
factorization framework [l3|, two-body Dirac equations of constraint dynamics [l4(], effec- 
tive Lagrangian (l5|, and bound state perturbation theory 16[. The lattice calculation 17| 



and rigorous QCD predictions 



191] were also applied. In addition, the s-wave and p- 



wave electromagnetic and light hadronic quarkonium decays in the heavy-quark velocity 
expansion were computed in the nonrelativistic QCD approach 20l-|22| . 

This paper is aimed at the study of the two-photon and two-gluon decay widths of p- 
wave heavy quarkonium states including the scalar (xco, Xho, x'm) ^^^ tensor (xc2, Xb2, x'hi. 
mesons. It is known that heavy quarkonium is relevant to nonrelativistic treatments [3 



231 ]. Although nonrelativistic QCD is a powerful theoretical tool used to separate high- 



energy modes from low-energy contributions, in most cases those attempting the calculation 
of low-energy hadronic matrix elements have relied on model-dependent nonperturbative 
methods. The light-front quark model (LFQM) 2J, |25|] is a relativistic quark model in 
which a consistent and fully relativistic treatment of quark spins and the center-of-mass 
motion can be carried out. This model has many advantages. For example, the light-front 
wave function is manifestly Lorentz invariant as it is expressed in terms of the momentum 
fraction variables in analog to the parton distributions in the infinite momentum frame. 
Moreover, hadron spin can also be correctly constructed using the so-called Melosh rotation 



261 ]. This model is very well suited for studying hadronic form factors. Specifically, as the 



recoil momentum increases (corresponding to a decreasing g^), we have to start seriously 
considering relativistic effects. In particular, at the maximum recoil point g^ = where 
the final-state particle could be highly relativistic, there is no reason to expect that the 



nonrelativistic quark model is still applicable. 



The LFQM has been employed to obtain some physical quantities 27H30|]. However, 
one often calculates a particular component (the plus component) of the associated current 
matrix element in the LFQM formulation. Because of the lack of relativistic covariance, the 
results may show some inconsistencies. The usual strategy of taking only the plus component 
of the current matrix elements will ignore the zero-mode contributions and render the matrix 
element noncovariant. As a consequence, it is desirable to construct a covariant light- front 
model that can provide a systematic way of exploring zero-mode effects. Such a covariant 
model has been constructed in [31i| for heavy mesons within the framework of heavy-quark 
effective theory. Without appealing to the heavy-quark limit, a covariant approach of the 



light-front model has been put forward for the usual s-wave mesons 32|, extended to the 
p-wave mesons [33| , and employed in the context of the s-wave heavy quarkonium 3J] . In 
this study, the p-wave heavy quarkonium is explored through this covariant model. The 
details and formalism are displayed in the next section. 

The remainder of this paper is organized as follows. In Sec. II, the formalism of a covariant 
light-front model is shown in cases of scalar and tensor quarkonia and the annihilation rates 
of these p-wave heavy quarkonium are derived. In Sec. Ill, after fixing the parameters 
which appear in the trial wave function, the numerical results and discussions are presented. 
Finally, conclusions are given in Sec. IV. 

II. FORMALISM OF A COVARIANT LIGHT-FRONT MODEL 
A. Formalism 

In the conventional light-front framework, the constituent quarks of the meson are re- 



quired to be on their mass shells (see Appendix A of Ref. 33| for an introduction). Jaus 



32| proposed a covariant light-front approach that permits a systematic way of dealing with 
zero-mode contributions. Physical quantities can be calculated in terms of Feynman momen- 
tum loop integrals which are manifestly covariant; this means that the constituent quarks 
of the bound state are off shell. In principle, this covariant approach will be useful if the 
vertex functions can be determined by solving the QCD bound state equation. In practice, 
we have to be content with phenomenological vertex functions, such as those employed in 



the conventional light-front model. Therefore, by using the light-front decomposition of the 
Feynman loop momentum, say p^, and integrating out the minus component of the loop 
momentum p~ , one moves from the covariant calculation to the light-front one. Moreover, 
the spectator quark is forced to be on its mass shell after p~ integration. Consequently, the 
covariant vertex functions can be replaced by the phenomenological light-front ones. 

As stated in passing, in going from the manifestly covariant Feynman integral to the 
light-front integral, the latter is no longer covariant as it can receive additional spurious 



contributions proportional to the lightlike four vector u^ 



{uj~,uj+,uj±) = (2,0,0±). The 



undesired spurious contributions can be eliminated by the inclusion of the zero-mode con- 
tribution which amounts to performing the p~ integration. The advantage of this covariant 
light-front framework is that it allows a systematic way of handling the zero-mode contri- 
butions and permits the user to obtain covariant matrix elements. 

To begin with, we consider the decay amplitudes given by one-loop diagrams, as shown 
in Fig. 1 for the two-photon decay of p- wave quarkonium states. The incoming meson has 
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FIG. 1: Feynman diagrams for the M — )■ 77* process where P in the parentheses denotes the 
momentum of the meson, (b) is related to (a) by the exchange of two photons. 



the momentum P = Pi + P2, where pi and p2 are the momenta of the off-shell quark and 
antiquark, respectively, with masses m. These momenta can be expressed in terms of the 
internal variables {xi,pj_), 



Pts = 3^1,2^"^, Pl,2± = Xi^2P± ± P±, 



(2.1) 



with xi -h X2 = 1. Note that we use P = (P+,p-,Pj_), where P^ = P'^ ± P^, so that 
P^ = P^P^ ~ Pi- III the covariant light-front approach, total four momentum is conserved 



at each vertex where quarks and antiquarks are off shelL However, it is useful to define some 
internal quantities analogous to those defined for on-shell quarks: 



"=+pi+P^, p. = (£lZ£l)«l, M„^ = (2e)^ = Pi+I^. (2.2) 

where m is the mass of heavy quark c or h. Here, Mq can be interpreted as the kinetic 
invariant mass squared of the incoming qq system and e, the energy of the quark i. 

We need Feynman rules for the meson-quark-antiquark vertices to calculate the ampli- 
tudes shown in Fig. 1. These Feynman rules for vertices of scalar (5*) and tensor (T) mesons 
are 

f -iHs, for 3Po, , , 

iTm = { . ^ (2.3) 

^t[7m - w-(Pi -P2)^{Pi -P2)v, for ^Pa- 



It has been shown in Ref. [35|] that one can pass to the light-front approach by integrating 
out the p^ component of the internal momentum in covariant Feynman momentum loop 
integrals. The specific form of the covariant vertex functions for on-shell quarks can be de- 



termined t 



irough a comparison with the conventional vertex functions as shown in Appendix 



A of Ref. 331] . They are obtained as 

Hm -^ Hm{P\,PI) = hM, 

Wv ^ Wv{pIpI) = wv, (2.4) 

where pi = P — P2 and P2 = m? ■ The form of the function Km contains two parts: one is 
the momentum distribution amplitude (t){xi,pi_) which is the central ingredient in light-front 
QCD, while the other is a spin wave function which constructs a state of definite spin (S*, Sz) 
out of light-front helicity eigenstates (Ai, A2). The spin wave function is constructed by using 



the Melosh transformation [26j and its spin structure is shown in Eq. (12. 3p . The explicit 
forms of Km and wy are given by 



IXiX2 1 Mq 



N^ V2Mo2y3 



1x1X2 1 



Wv = Mo + 2m. (2.5) 

The momentum distribution amplitude 4>{xi,p±) is the generalization of the distribution 
amplitude 4>{xi) of the pQCD method and can be chosen to be normalizable, i.e., it satisfies 



dxd p± 
2(27r)3 



x,p^)\' = l. (2.6) 



In principle, (f){x,p±) is obtained by solving the light-front QCD bound state equation 
HlfI"^) = M|\l/) which is the familiar Schrodinger equation in ordinary quantum mechanics 
and Hlp is the light-front Hamiltonian. However, at the present time, methods concerning 
how one can solve bound states equations are still unknown. We are satisfied with utiliz- 
ing some phenomenological momentum distribution amplitudes which have been constructed 
phenomenologically in describing hadrons. One widely used form is the Gaussian type which 
we will employ in the application of the covariant light-front approach. 

B. Two-photon decay widths of p-v/ave heavy quarkonium 

Quarkonia are eigenstates of the charge conjugation operator C with eigenvalues C = 
(— l)-^^'^, and the charge conservation requires charge conjugation C = +1 state coupling 
to two real photons. Thus, only ^Pq and ^P2 levels of the p-wave quarkonium states can 
transform into two real photons. In the process, M — )• 77 , the final two photons are both 
on shell. For the purpose of illustration, it is useful to consider a more general process: 
M — )■ 77* with one photon off shell. We introduce some transition form factors fi{q^) 
arising from the M77* vertex. The M — )■ 77 process is related to the form factors at q^ = 0, 
i.e., /j(0). It is worth mentioning that, however, the state M in the process M — )• 77* is not 
only the ^Pq and ^P2 levels. 

The matrix elements of the decay S — )■ 77* and T — )■ 77* have the following structures 

u 



Af, ^ {j{k,e)\qj,q\S{P)) = fo{q')[e,{P-k)-k,{P.e)], (2.7) 

Al = (7(A;,e)|g7^g|r(P,e)) 

= fiil'^) i'^M -e-k)- k^{q ■e-e)+ e^^e^iq ■ k) - e^^k^ie ■ q)] 

+ f2{q')h{q-k)-k,{e-q)]^^, (2.8) 

q ' K 

where e is the polarization vector of the on-shell photon, q = P — k is the momentum 
transfer, and e is the polarization of the tensor meson which satisfies the following relations: 

P%,(A) = P%,(A) = 0, £^.(A) = £,^(A), £^^(A) = 0, £^^(A)£;,(A')=<5aa', 
Y, £^,(A)4^(A) = -M^aM.p + -M^pM,^ - -M^^M^p, (2.9) 

A 



with A = ±2, ±1,0 representing the tensor meson hehcities and M^y = g^^ — P^P^/M"^. 
There is no exphcit representation of the meson polarization tensor £:^;y(A). The transition 
amphtude for the processes of S* — )■ 77* and T -^ 77* can be derived from the common 
Feynman rules and the vertices for the meson-quark- antiquark coupling given in Eq. f l2.3p . 
In the covariant light-front approach, the meson is on shell while the constituent quarks 
are off shell. To the lowest order approximation, S(T) — )■ 77* is a one-loop diagram and is 
depicted in Fig. 1. The amplitude is given as a momentum integral 

+ ^^^^^J r [(- ^2 + m)7^(;^6 + m) /(:^i + m)] |, (2.10) 
Al = ^eyN^^^ I ^^^^^Tr [t{- ^2 + m) /(^, + m)7^(^i + m)] 



where 



+ ^^^^^J r [t(- ^2 + m)7^(^b + m) /{^^ +m)]\, (2.11) 



Q-P2, 

2 2 I • 

= P2 — m + te, 

2 2 I • 

= p^ -m +te, 



Pa = Pi- q, 


Pb-- 


Ni= p\ — ni^ + ie. 


N2 


Na = pl-m^ + ie, 


Nb 



t = e 



a/3 



r-^{pi-P2r 



{Pl - P2 



1/3 



(2.12) 



and Cg is the electric charge of the quark: e^ = 2/3 for the c quark and Cg = —1/3 for the b 
quark. The first and second terms in Eq. ( ]2.10p come from diagrams Figs. 1 (a) and 1 (b), 
respectively. In the calculation, it is convenient to choose the purely transverse frame g"^ = 0, 
i.e., g^ = —q'j_ < 0. The advantage of this choice is that there is no so-called Z- diagram 
contributions. The sacrifice associated with this approach is that only the form factor at 
spacelike regions can be calculated directly. The values at the timelike momentum transfer 
g^ > regions are obtained through analytic continuation. In this study, the continuation 
is not necessary because we only need the form factors at g^ = for the S{T) — )■ 77 process. 
First, we discuss the calculation of Fig. 1 (a). The factors A^i, A'2, and Na produce three 
singularities in the pj~ complex plane: one lies in the upper plane and the other two are in the 
lower plane. By closing the contour in the upper pi complex plane, the momentum integral 
can be easily calculated since there is only one singularity in the plane. This corresponds 



to putting the antiquark on the mass shelL Given this restriction, the momentum p2 -^ P2 
with P2 = m^ 1 and pi = P — P2- The on-shell restriction and the requirement of covariance 
lead to the following replacements: 



iVi 



iVi =xi{M^ 



M, 



07' 



N2 ^ N2 = Ni + {l- 2xi)M^ = X2M^ - xiMl, 
Na ^ Na = X2q^ - XiMl + 2px ■ g±, 



(Tpi Hs,T 



-)■ 



-?7r 



dx2d p± hs^T 



(2.13) 



{2Tl*)N^N2Na J (27r4) X2NiNa 

For Fig. 1 (b), the contour is closed in the lower pi complex plane. It corresponds to putting 
the quark on the mass shell and the momentum pi — )■ pi with p\ = m^. For this scenario, 
we need to do the following replacements: 



N^ ^ Ni= xiM^ - X2M^xi{M^ 

N2 -^ N2=X2{M^ 



Ml). 



Ml), 



Nb ^ Na = xiq^ - X2MI - 2px ■ g±, 



d'^Pi Hs,T 



-)■ 



-tTT 



dx2d p± hs T 



(2.14) 



(27r4) iViiVsiVb ' ^" J (27r4) xiiYsiY,' 
After the above treatments, the transition amplitudes of 5* — ?■ 77* and T — )■ 77* for Fig. 1 
(a), for example, are obtained as 

dx2d'^p± hs m Mq 



^S{a) 



N.ey 



An^ XiX2(M2 



M2 



-X2g2 + xiMl - 2p^_ ■ q^ 2^3 



X {e^{m^ -P-q-pl + 2pi-q) + 2pi^{2e ■ pi - t ■ P - e- q) + q^{e ■ P - 2e ■ pi] 



+P^e ■ g}. 



jT(a) 






(2.15) 



dx2d'^pj 
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47r3 xiX2iM^ - Ml) -X2q^ + XiMl - 2p± ■ q± 
m + P ■ pi — Pi)pi ■ e ■ q + [m — P ■ q + 2pi ■ q — p^) ( 1 



2m 

Wv 



Pi-e-Pi 



+P, 



[m 



Pi- q- Pi)e ■€ -pi- e- pipi ■€ ■ q + e- qil 



2m 



Wv) 



Pi-e-pi 



{p\ — m? — P ■ pi)t ■ e ■ pi + {t ■ P — 2t ■ pi) ( 1 



2m, 

Wv 



Pi-e-pi 



[m 



P ■ q + 2P ■ pi — p\)e ■ e ■ pi + e ■ Ppi ■ e ■ q 



-(4e-pi-2e-g-2e-P) (l 



2m\ 

Wv) 



Pi-e -pi 



+£f,uPi 



m, {e ■ pi — e ■ P — e ■ q) + {e ■ piP ■ q — e ■ Ppi ■ q — e ■ qP ■ pi) 



+Pi(e ■ q + e- P - e- pi) 



(2.16) 



The integration of pi^, pi/^piu , Pi^^PiuPia, and Pi^PiuPiaPip in Eqs. (I2.15P and (I2.16P can be 
expressed in terms of three external vectors: P, g, and w, as in the Appendix. The transition 
amphtude of S{T) — )■ 77* for Fig. 1 (b) can be obtained through a similar process. If we 
choose the frame where the meson is at rest and the photons travel in the ±2; directions, 
then the two-photon decay amplitude of S" — )■ 77 is obtained as 



A^^(S ^ 77) = (^f^+^f^) 1.^=0 -e' 



n o Nc , f dx2d p± mc 
e„e \ — e-e / ; 



' 2xiM^ - M2 - 4pl 2x2M^ -M^- 4pl \ 



m'^ + Pj_ 



2^/?>XiMq 



2^X2Mq 



where e' is the polarization of another photon. The decay rate of the process S* — )■ 77 is 36 1 



T{S ^ 77) 



8^^l pel 



k\Y.\M'\\ 



(2.18) 



— * 

where s = 1/2 for two identical photons and the photon momentum \k\ = M/2. Regarding 
the case of T — )■ 77, in the practical calculations below, it is convenient to represent e^^ in 
the following form 37l |: 



£^,u{>^) 



v/6 
1 



(2-|A|)(l-|A|) 



3n^n^ +[g^,u- -j^ 



+-(1 - |A|) [A(« - n^O + ^|A|(« + nX)] 
+1(2 - |A|) [A(n^^ + nX) + ^|A|(n^< + n^)] , 



(2.19) 



where n^{i = 1,2, 3) is the basis vector satisfying n^n'^g^i, = Qij. In the center-of-mass frame, 
we may have P'^ = (£',0,0, |P|) and choose the basis with coUinear n-^ and P^ vectors as 
the simplest one 

1 



< = (0,1, 0,0), n^ = (0,0,1,0), n^ 



M 



(|P|,0,0,E). 



(2.20) 



For the frame in which the meson is at rest and the photons travel in the ±z directions, the 
decay amplitudes of T — )■ 77 are obtained for A = and ±2 
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(2.21) 
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^^ Jm^ + pI \XiM^ X2M( 
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(2.22) 



where 



t'Oa 



xi(m2 - X2M^){xiMl - X2M^) + p\[-?,m^ + (1 - 4xi + 2xJ)M^ 



P± 
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+ — {xi{2m^ - M^){x2M^ - xiMl) + p\[2m^ - ixlM^ + (4xiX2 - l)M^] - 2pi} 

Wv 



t, 



06 



I 



Oa 1x1-0X2) 



m 



t±2a = p\\m^ + MqX? + 2 p] I (ze' ■ ^26 -ni^-ie- n2e' ■ rii ± e ■ nie ■ ni ^ e ■ n2e ■ 77,2), 

V Wv J 

t±2b ^ t2a\x-i-^X2- 

The derivations of Eqs. fl2.2ip and (12.221) use the formulas 



2.23) 
(2.24) 



0. 



P±9± n 3p±g± 



,0, 



. (2.25) 



j d^p± \pi. ■ q±, (p± ■ g±)^ (p± ■ g±)^ (p± ■ gx)"^] = j d'^p^ 

The amphtudes of T — )■ 77 for A = ±1 vanish because the combination of the hehcities of 
two final photons can be ±2 or 0, but never equal to ±1. The decay rate of the process 

T — )• 77 is 

S \k\ y;^ V^ I A aT\2 



r(T -^ 77) 



E j:\m' 



SvrMi 5 ;,=o,±i,±2 pel 



(2.26) 



The factor 5 in the denominator corresponds to 2J + 1, where J is the total angular mo- 
mentum of the meson. 

Finally, the two-gluon decay width of quarkonium can be easily obtained from the two- 
photon decay width, with a simple replacement in the photon decay width formula 

2 
— ( 
9 



e^a^ — >■ -at- 



(2.27) 



III. NUMERICAL RESULTS AND DISCUSSIONS 



In this section, the two-photon and two-gluon decay widths of p-wave heavy quarkonium 
states are estimated. Prior to numerical calculations, the parameters rric^b and Pccbby which 
appeared in the wave function, must be first determined. We consider the Hamiltonian of 
the s- and p-wave heavy quarkonium states as 



4a „ 



Hs = 2^/rn? + n? + Konf - -77^ + QoSi ■ S2, 

Sr 

4a 

Hp = 2Vm2 + K^ + Vconf - ^ + giS ■ L + g2Si2 + gsSi ■ 

3r 

10 



S2, 



(3.1) 
(3.2) 



where Konf = 6r(6r2) is the hnear [harmonic oscillator(HO)] potential, S'12 = (3si -f S2-r — 
si ■ S2) is the tensor force operator, and (70,1,2,3 are the functions of the relevant interquark 



potentials (the details are shown in, for example, [38|, l39|]). In this way, not only is the 
spin-weighted average of the s-wave states M{Sj) = [M{^Sq) + 3M(^5'i)]/4 free of the spin- 
spin contribution, but also the mass difference between the spin-single ground state M(^Pi) 
and the spin-weighted average of the p-wave triplet states M(^Pj) = [M(^Pq) + 3M{^Pi) + 
5M{^P2)]/9 has a contribution which comes from the spin-spin interaction.^ Experimentally 
the latter hyperfine splitting is less than 1 Mev [1| in the charmonium sector and can be 
ignored here. Notably, we can use the masses M{Sj), Mi^Pj) and their variational principle 
for the Hamiltonian equations. (13.1 p and (13. 2p to determine parameters m and /3. In the 
process, the 15, IP, and 2P states harmonic wave functions 

3/4 



(a;,P±) 
0^^(a;,p±) 



A\^ 






4 Ul, 



/3^ 



4 



7r 



3/4 



3/4 
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2/32 _ 


7 






dpz [Y 


2/32 _ 


; 




dp. [^ (^ 
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5f3\ 
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2/32 _ 



(3.3) 
(3.4) 
(3.5) 



which satisfy the normalization Eq. (12. 6 p and their conjugate coordinate wave function 

52 \ 3/4 



a^/ 



71 



exp 



2^2 



(3'r 



^i:{r) 



k2Pl 



/^2\3/4 

V2 ( — 1 /3^m exp 



2^2 



/3V 
2/32r2^ 



exp 



2^2 



P'r 



(3.6) 
(3.7) 
(3.8) 



where am=±i = T{cix ± "i^y)!^^ and 0^=0 = dz are needed. In addition to the coefficient of 
confined potential 6, the cc and 66 sectors each have four parameters, tti^, a^, /3^'^, and /3^^, for 
the 15* and IP quarkonium states. Regarding the constraints, in addition to the four masses 
M^^i^liXSj) and iV4gf,j;(l^Pj), the four equations (13. 3p . (13. 6p and (13.40 . (13. 7p are used as the 
trial functions of the variational principle for IS* and IP states, respectively. Therefore, h is 
the only free parameter in our ffiting. We employ the data r(xco ~^ 77) = 2.36 ± 0.35 keV 



^ The calculations of the expectation values of the fourth and fifth terms for the Hamiltonian Eq. p.2p can 
be referred to in the appendix of Ref. [40| . 
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TABLE I: The relevant parameters b, nig, a^, and /3 of the p-wave heavy quarkonium states. 



Potential 


b 




rUg (GeV) 


as 


Z?!^ (GeV) 


/32-P( GeV) 


Linear 


O.mt'oZ GeV2 


cc 
bb 


L42^0.02 
4.78^0.01 


0.489;0:0}« 
n -^00-0006 

u.oyy_,_0.007 


0.510 ±0.002 
0.807 ± 0.005 


0.568_o:o24 


HO 


O.Omt^,Zl GeV3 


cc 

bb 


1.42:^0.02 

A 07-0.01 

4-«'+o.oo 


n Qt;o-0.025 
U'-JJO+O-OIS 

r, .nr:-0.005 

U.y:ZO^QQQg 


n C77+0.001 

'-'•"-" '-0.003 

0.847 ± 0.05 


n c,7-| +0.007 
'^•^'^-0.009 



fl 



Jj, to fix 6; then the above parameters can be determined and shown in Table 1. For 2P 
states, due to the insufficient data regarding x'cj^ o^ly the parameter /3^^^ is determined by 
the mass Mi^i{2^Pj) and also revealed in Table 1. There are three items in Table 1 worth 
mentioning: First, the parameter b = 0.176io;oo8 GeV^ in the linear potential is consistent 
with the string tension 6 = 0.18 GeV^ which is well known from other quark model analyses 
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4l|,U- 



Second, there are presently different conceptions about the value of a^ in the 



low-energy region. In lattice QCD 



and the field correlator method 



for example, ones 



found for the coupling constant in the static potential, parametrized as a linear plus Coulomb 
potential, the small values a^ are 0.21 and 0.16, respectively. In phenomenological potentials 
the Coulomb constant is larger: ag is 0.46 4^ and 0.43 ± 0.02 451] which corresponded to 



the charmonium and bottomonium states. Here our values are consistent with the latter. 
Third, in the HO potential, the strong coupling constant a^ of the cc sector is smaller than 
that of the bb sector, which violates the concept of asymptotic freedom. The dependences of 
as{cc) and as{bb) on b regarding the linear and HO potentials are shown in Figs. 2 and 3, 
respectively. Obviously, the HO potential is not suitable when b > 0.0410 GeV'^. However, 
in the data fitting of r(xco ~^ 77); we obtained b = 0.0490;l^o;oo24 GeV'^. Therefore, only the 
linear potential is employed in the following calculations. 

Next, we use the parameters of the linear potential in Table I and Eqs. fl2.18p . (12.261) . 
and (I2.27P to calculate the two-photon and two-gluon decay widths of the Xco,c2, Xbo,b2, 
and x'bo 62 states. The numerical results, which compare with the experimental data and 
other theoretical evaluations, are shown in Tables 2 and 3. From these tables, except in 
the case of our two-photon decay width of Xco as an input, all the cc decay widths of this 
work are consistent with those of previous experiments and the major theoretical methods. 
However, our results of Xbo,b2 and Xw 62 differ from other calculations. More specifically, our 
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FIG. 2: Dependences of as{cc) (solid line) and as{bb) (dashed line) on b in the linear potential 
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FIG. 3: Dependences of as{cc) (solid line) and as{bb) (dashed line) on b in the HO potential. 

two-photon decay widths of Xbo,b2 and Xbob2 ^^^ significantly smaller than those of other 
approaches. From Table 2, the ratio T{xco ~^ 77)/r(xfeo ~^ 77) is about 300 for our work 
and about 60 ~ 80 for other estimations. It is well known that the decay rate has two 
contributions: one is the kinematic phase space; the other is the dynamic decay amplitude 
square. No matter what approach, the phase space and factor e^ of T{xco -^ ll)^ in total, 
is about 50 times the one of r(xf,o -^ 77)- Regarding the comparison of decay amplitude 
}A^ in Eq. (I2.17p . we find that the contribution of the numerator in the parentheses is 
roughly the same for Xco and Xbo- On the other hand, the denominator in the parentheses. 
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TABLE II: Two-photon widths of the p-wave heavy quarkonium states. R is defined as r2^^/r2^°. 





r^- (keV) 


r^- (keV) 


R 


r^- (eV) 


r^- (eV) 


r^;" (eV) 


^i' (eV) 


PDG [1] 


2.36 ±0.35 


0.515 ±0.060 


0.218 ±0.058 










This work 


2.36 ±0.35 


0.346iS°?? 


0.147 ±0.026 


8.02i-^ 


3.99±riJ 


2.26tri? 


1 79+0.39 
^•'^-0.33 


Miinz [5] 


1.39 ±0.16 


0.44 ±0.14 


^9+016 


24 ±3 


5.6 ±0.6 


26 ±2 


6.8 ±1.0 


Wang [6 , 7] 


3.78 


0.501 


0.133 


48.8 


7.4 


50.3 


7.7 


Ebert [10] 


2.9 


0.50 


0.17 


38 


8 


29 


6 


Laverty'' [8] 


2.02(2.12) 


0.46(0.19) 


0.23(0.09) 


32.9(94.4) 


7.19(5.38) 


34.1(94.5) 


7.59(5.57) 


Schuler [13] 


2.5 


0.28 


0.11 


43 


7.4 






Gupta'' [9] 


6.38(8.13) 


0.57(1.14) 


0.09(0.14) 


80(85) 


8(12) 






Huang [18] 


3.72 ±1.11 


0.49 ±0.15 


n+oii 










Bodwin [19] 


7.1 ±2.5 


0.81 ±0.29 


0.11 ±0.08 










Crater'^ [14] 


3.96(3.34) 


0.743(0.435) 


0.188(0.130) 










Barbieri [4] 


3.5 


0.93 


0.27 










Godfrey [^2] 


1.29 


0.46 


0.36 










Lansberg [15] 


5.00 


0.70 


0.14 










Dudek*^ [17] 


2.41 ±0.58 














Lakhina [16] 


3.28 














LO(NLO)'= [46] 






0.267(0.184) 











"The values are obtained by the perturbative (nonperturbative) calculation. 

''The values are obtained by the QCD potential (alternative treatment). 

"^The values are obtained by the two-body (naive) decay amplitude. 

''Only the statistical error is shown. 

'^These values refer to a leading (next-to-leading) order calculation done at the renormalization scale 2mc- 

Mq, is approximately equal to the meson mass for the heavy quarkonium. Therefore, the 
dynamic decay amplitude square of Xco — ^ 77 is about 6 times that of xto -^ 77- Our ratio 
r(Xco -^ 77)/r(x6o -^ 77) — 300 is obtained by combining the above two components. 
Regardless, more experimental measurements are required for these channels. Finally, the 
ratio R = Fg^^^/Fg.^'' of experiment and various theoretical estimations are also listed in Table 
2 for a comparison. 
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TABLE III: Two-gluon decay widths of the p-wave heavy quarkonium states. 





rf;° (MeV) 


Tlf (MeV) 


rg° (keV) 


Tf^ (keV) 


rgo (keV) 


rg^ (keV) 


PDG" [IJ 


10.4 ±0.7 


1.98 ±0.11 










This work 


11 q+o-^ 


1 74-0.08 


43111E 


214;? 


1221^ 


92.3liL^ 


Wang[6,:a 


10.3 


2.64 


887 


220 


914 


248 


Laverty'' [8] 


4.68(4.88) 


1.72(0.69) 


960(2740) 


330(250) 


990(2740) 


350(260) 


Gupta'' [9] 


13.44(17.10) 


1.20(2.39) 


2150(2290) 


220(330) 






Bodwin [19] 


4.8 ±0.7 


1.98 ±0.18 










Barbieri [4] 


2.4 


0.64 










Godfrey [12j 


6.25 


0.774 


672 


123 


672 


137 


Ebert [llj 






653 


109 


431 


76 



i tot = i 2g • 

^The values are obtained by the perturbative (nonperturbative) calculation. 
"^The values are obtained by the QCD potential (alternative treatment). 



IV. CONCLUSIONS 

This study has discussed two-photon and two-gluon decay widths of p-wave heavy quarko- 
nium states within the covariant hght-front quark model. This formalism, which preserves 
the Lorentz covariance in the light-front framework, was applied to annihilations of the scalar 
and tensor quarkonia. To obtain the numerical results, we used the harmonic wave func- 
tions and fixed the parameters appearing in them. The constraints were the spin-weighted 
average masses M{Sj, ^Pj) and their variational principle regarding the Hamiltonian. We 
considered the linear and HO potentials in the Hamiltonian and found that, when the data 
r(Xco -^ ll) were fitted, the former resulted in a value of h consistent with that of other 
quark models and the latter led to a violation of asymptotic freedom. Therefore, only the 
parameters corresponding to the linear potential were applied to estimate the relevant decay 
widths. The numerical results showed that, for the cc sector, all of the decay widths were 
in agreement with the experimental data and the major theoretical calculations. However, 
for the bb sector, discrepancies appeared in the decay widths of Xbo,b2 and Xm b2 fi'om other 
estimations. 
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Appendix A: Formulas for the product of several pi's 

In general, after the pi integration, pi can be expressed in terms of three external vectors, 
P, q, and u, where P = P + k. Furthermore, the inclusion of the zero- mode contribution 
cancels the u dependence and in practice for pi in the trace under the integration, we have 



32 



i(i)^^ aW 



p,, = P,AY> + q,A',", ^ ^ (Al) 

Pi,Piu = g,uA?^ + P.PuA?^ + ih^u + q,Pu)Af^ + g.g.Af + ^^^'^+>^'^ i?jfA2) 

uP 

PlliPluPla = iOi^uPa + 9fiaPu + 9uaPfi)Af^ + {g^uQa + Qi^aQu + guaqii)A^2^ + P^,P^PaAf^ 



HP^^P.qa + P^QuPa + q^.PuPc.)Af + {q^q.P^ + g^P,g« + P^q,q^)Af 
+qi.quqaAf + -^{P^P.Ua + P^^UJ^Pa + UJ^P,P^)Bf'^ 

uP 

H ~[{PfJ.q<^ + qfj.Pu)(^a + {Pf^qa + qfj.Pa)uJu{Puqa + q^Pa)(^fj]Bi^\ (A3) 

uP 

9 4 

Pl^PliyplaPllB = 51 ^^M^^"/? A^"^^ + H ^iM!^"/3^f ^ (^4) 

where the symbol = denotes that the equation is true only in the integration and 

hfiuap = g/iugap + g^iagvp + gti^g^a, 

hfiuafB = gnuPaPp + gfiaPi/P/S + g/i/sPuPa + gaf^P/iPu + gi^/^PfiPa + guaPfiPpi 

hfMuaf} = gfiuiP^qp + Ppqa) + permutations, 

hfiuap = g/iuqaqp + g^iaqvqp + g^ipquqa + gapq/iqu + gupq^iqa + guaq^iqp, 

hfiuap = P/iPiyPaqfi + PfiPuqaPfB + PfiquPaPp + qfiPiyPaPp, 

hfiuap = PfiPiyqaqiB + pcrmutatious, 

hfiuajs = Piiqiyqaqi3 + q^Puqaq^ + qfiqyPaqp + q^quqaPpi 
ht^uaii = qfiquqaqpi 

Jiixvap = -^[gf^L'iPai^f} + Pp^a) + permutatious] , 
uP 
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with 



(1 



A 



(2 



A 



(3 



(3 



A 



Af 
4^' 



B. 



(4: 






Xi 

y 



— -^iPfj.PyPa^ 13 + PfiPu^aPli + Pfi^uPaP/S + ^fiPuPaP/s), 

uP 

1 



^[(-P/x-Pi^ga + P^j,quPa + q^,PyPa)(^p + permutations], 

^[(-PM^i^iJa + ^M-Pi^^a + QtiQuPa)^!^ + permutations], 
W-P 



(A5) 



A« 



.(1) P±-9± .(2) 

^1 ^5 ' ^1 



2 (P± ■ ?±> 



g^ 



Af 



(^i 



{1)n2 



^(1)4(1) 4 



(2) 



{A 



(1)^2 ^1 



(2) 



g^ 



B\ 



(2) 



4^^Z2-4'\ 4'^ = 4^M 



1 ) 



l(l)/l(2) 



I (3) 



(1) /l(2) 



4^ M^ -" 4^'''' — A^ > A^ > A^^^> — A^ ' A 

^2 ^1 ; ^3 ~ ^1 ^2 ; ^4 — ^2 ^ 



l(3) 



l(l)/l(2) 



2 ; 



A 



(3) 



4(1)^(2) 



41)42) _ ^^(1)^(2)^ 5(3) ^ ^(1)(5(2) _ ^(2)^ 



5(^) 



Ai^Bf^ + ^-^Af\ 



q- 



{A\ 



(2)n2 



A 



(4) 



4(1)4(3) 4(4) _ 4(1)4(3) 4(4) _ 4(1)4(3) Ai 

/il n.1 , /I3 — n.1 /±2 ; ^4 ~ ^2 ^2 ~ 



(4) 



gZ 



4(1)4(3) 4(4) 

^1 ^3 ) ^6 



(1)4(3) 



Ai^M, 



2.4(4) 
,2^4 5 



4(1)4(3) 

?(4) 



4^^ = A\'>A 
(1)4(2)^ 



(1)4(3) 4(4) 

5 5 ^8 



4(1)4(3) 
^1 ^6 5 



5r = Ai^Mi^^z2-Ai 



1(4) 



B. 



(4) 



Al^^i?}^^ - A 



2 ! 



A(^)i?f - 4^ 5f = ^!'^ f4^^^2 + 2^-^A^;^A?A - 4'), 



Z2 = iVi + (X2 - Xi)M2 + (g2 + ^ . P) 



g2 



(A6) 



Equation flA6p is obtained by contracting the larger number of pi's with u^, qis, and g^/s, 
and through a comparison with the complete expression of the fewer pi's. 
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